Global convective cell formation in pair-ion plasmas 
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Abstract: The global electrostatic mode in pair-ion plasmas is discussed for cylindric 
geometry and for a radially inhomogeneous equilibrium density. In the case of a Gaussian 
radial density profile, exact analytical eigen solutions are found in terms of the Kummer 
confluent hypergeometric functions. The mode is identified as a convective cell propagat- 
ing in the poloidal and axial direction, having at the same time a standing wave structure 
in the radial direction. 



PACS Numbers: 52.27.Cm; 52.30.Ex; 52.35.Fp 



1 



A new research field and a new area of increasing scientific activity has emerged re- 
cently after a series of experiments^"^ in which a pure pair-ion plasma has been produced, 
i.e., a plasma without electrons. The two ions, C^, are produced in a simultaneous pro- 
cess of impact ionization and electron attachment, and further collected by a magnetic 
filtering effect, i.e., by a diffusion in the radial direction (perpendicular to the magnetic 
field lines). Being separated from electrons, the ions are then collected in a narrow and 
elongated chamber (90 cm long and with a diameter of 8 cm). In Refs. 2-4 several types 
of modes have been reported in such plasmas, viz. the ion plasma wave (IPW), the ion 
acoustic wave (lAW), and the so-called intermediate frequency (IF) wave. The detailed 
measurements presented in the most recent Refs. 3 and 4, reveal that the I AW actually 
has two separate branches, and is accompanied by some additional backward propagating 
mode between these two branches. Some features of the experimentally observed modes 
still remain unexplained. 

These experimental works have been followed by numerous studies^"^^, dealing with 
various aspects of linear and nonlinear waves and instabilities in pair-ion plasmas. In 
addition to this, the recent successful production of a pair-hydrogen plasma,'* that is even 
more attractive because of the small ion mass and the obvious consequences related to 
this, indicates that the number of experiments and analytical studies in this particular 
field will grow even further. The topic is interesting also in view of the fact that a similar 
pair plasma, comprising much lighter particles (electrons and positrons) in the past years 
has also been created under laboratory conditions. The knowledge of processes in such 
plasmas will help us to understand the physics of some pair-plasmas in space, comprising 
electrons and positrons, like in the atmospheres of pulsars and in active galactic nuclei, 
and even in flares in the lower solar atmosphere. 

So far, studies dealing with waves and instabilities^^^* have been carried out mainly 
in the framework of a local mode analysis and in Cartesian geometry. However, the 
experimental conditions mentioned above^"^ are such that, in some domains of frequencies 
and wave lengths, the proper cylindric geometry is to be used and the effects of the plasma 
boundary on the modes are to be taken into account. Therefore, in the present work we 
shall focus on frequencies below the ion gyro-frequency, on the effect of radial density 
variation (and the corresponding radial dependence of the mode amplitude), and on the 
effect of the boundary on the spatial structure of the mode. 

We consider a cold cylindric pair-ion plasma configuration with an externally given 
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magnetic field in the axial direction Bq = BqCz, and with a radially varying density 
profile nQ{r). We note that the spatial variation of the density in both the radial and 
axial directions has been experimentally detected/ with the radial density profile being 
nearly Gaussian. Such equilibrium density gradients are in fact a rather common feature 
of various laboratory plasmas. Note that an axial density inhomogeneity has also 
been detected. Yet, in view of very different axial and radial lengths of the chamber, that 
inhomogeneity may be neglected without losing much of essential physics. 

In the presence of electrostatic perturbations of the form ~ /(r) exp{—iujt+ik-f+ikzz), 
propagating at an angle with respect to the magnetic field vector, from the momentum 
equations for the two species (j = a, b) we obtain 
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Here, Q = eBQ/m, rua = rrib = m, the — sign stands for positively charged ions, and the 
other notation is standard. Equation ([T]) is used in the two continuity equations which, 
with the help of the quasi-neutrality condition, yield 

For perturbations having a wave number component in the poloidal direction, we have 
k ■ r ^ mO and from Eq. ([2]) we obtain a differential equation for the amplitude of the 
perturbed potential 



/I 1 rfr2o\ d ^2 i^^ ^ 



0. (3) 



Without the density gradient, the general solution of Eq. ([3]) can be written in terms 
of the Bessel functions of the first kind, i.e., 0(r) = ciJm[b^^'^r] + C2J-m[b^^'^r], where 
b = k1{VL^ /uj"^ — 1) denotes the eigenvalue. For a plasma extending in the radial direction 
up to r = i?, the proper boundary condition requires vanishing solutions at the boundary, 
so that fe^/^i? = ei gives the dispersion equation. Here, e; is the /-th zero of the Bessel 
function J±m- 

In the presence of the density gradient, and assuming a Gaussian density profile 
no(r) = Ai'o exp(— K^r^), the eigenmode equation becomes 



4>{r) = 0. (4) 
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The general solution of Eq. (jl]) is^^'^^ 

b m 



b m 2 2 

1 — m, K^r 



■ (5) 



Here, iFi is the Kummer confluent hypergeometric function, and the Ci^2 correspond to 
the integration constants. 

The eigenvalue b yields the eigenfrequency 

which describes the obliquely propagating (in the 9 and z directions) electrostatic con- 
vective cells, i.e., modes twisted around the magnetic field vector. Examples of similar 
convective cells in the domain of the ion gyro-frequency, and the lower-hybrid frequency, 
are known from the literature. 

Remarkable properties of the pair-ion plasma can be seen by comparing Eqs. and 
in]) with the corresponding equations in electron-ion plasmas, later describing the drift 
wave. In the present case, the poloidal wave number m directly determines only the 
radial eigenfunctions. The same holds for the equilibrium density profile parameter k. 
The dispersive properties of the mode are directly determined only by the eigenvalues b. 
However, in the electron-ion plasma^^ the dispersion equation contains all the parameters, 
viz. the poloidal wave number m, the density parameter k, and the eigenvalue b. The 
difference between the two cases appears due to the pair property of the plasma in the 
present case, which includes the cancelation of terms with the E x 5-drift, i.e., the terms 
veVho, in two combined continuity equations. 

Well behaving solutions are finite in the plasma column and should not have a finite 
radial velocity aX r = R and at r = 0. Hence, d(j)i{r = 0,l)/d6 = 0, and this is satisfied 
with a potential profile vanishing at the axis and the boundary. Therefore, we set Ci = 0, 
and the appropriate profile will be found for certain values of the parameters k and b. 

Observe that the function iFi[di,d2,r] is a polynomial with a finite number of terms 
if di < 0, and if in the same time (^2 > or (i2 < c^i. This implies oscillatory (standing 
wave) solutions in the radial direction. In the present case, we have d2 = 1 + m > 0, and 
oscillatory solutions di < exist if 6 > 2mK,'^. For example, iFi[-l,2,r2] = 1 -rV2, 
iFi [-2, 2, r^] = 1 - r2 + iFi [-3, 2, r^] = 1 - 3r^/2 + - rV24. 

The eigen- values b satisfying the boundary conditions (e.g., a vanishing potential at 
r = R) can be found numerically. There exists a multiple choice of eigenvalues 6, in terms 



of K, as shown in our recent Refs. 21, 22. The values of k may be taken from the interval 

0. 1 to 1.5. It is easily seen that for the assumed Gaussian density n{K,, r) normalized to Nq, 
where r is normalized to R, this yields ?t,(0.1, 1) = 0.99 and n{1.5, 1) = 0.1, respectively. 
For this range of k, and taking as example the poloidal mode number m = 3, the first 
branch of the eigenvalues b in terms of k is presented in Fig. 1. These pairs of b, k, give 
the first eigenfunction ■ iFi, i.e., the profile with the first node of the radially oscillating 
standing wave solution ai r = R. It can be shown that for m = 3 and for k, = 0.1, the 
second branch starts with b = 95.3. This branch gives the second eigenfunction, the one 
with the second node of ■ iFi at r = R. The third branch begins with b = 169.2, etc. 

The three first eigenf unctions (standing wave solutions in the radial direction) are 
presented in Fig. 2 for k = 0.3, and b = 40.53, 95.18, 169.2, respectively. A full three 
dimensional plot of the modes gives twisted solutions of the form cos{m6+kzZ—u!t)-r"^-iFi, 

1. e., electrostatic waves traveling in the poloidal and axial directions, having at the same 
time a standing wave structure in the radial direction. 

The dispersion equation, Eq. shows that the frequency of the mode remains 

below the ion gyro-frequency. Each discrete eigenvalue b gives a separate dispersion line, 
the first three lines being presented in Fig. 3. Hence, for the given density profile with 
K = 0.3, any of these branches of oscillations can take place. 

To conclude, in this Brief Communication exact analytical solutions are presented, 
that should be expected in the low frequency range below the ion gyro-frequency in the 
recently produced radially inhomogeneous pair-ion plasma. Contrary to the similar case 
in electron-ion plasma, where in this frequency range the drift mode is found, ^^'^^ in the 
present case, because of the pair properties of the plasma components, these modes are 
convective cells, and the mode frequency is not directly determined by the density gradi- 
ent. Global modes determine global properties of a plasma, so that after the identification 
of the modes described here, their eventual observation should be used in the diagnostics 
of pair-ion plasmas. 
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Figure captions: 

Fig. 1. Locus of pairs b, k for tlie first eigenf unction, satisfying tlie condition ■ iFi = 
at r = i?. 

Fig. 2. The three lowest radial eigenfunctions • iFi, in arbitrary units, for k = 0.3. 

Fig. 3. Normahzed frequency ([6]) for the three first eigenfunctions from Fig. 2, for the 
given density profile with k = 0.3. 
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